Adhesion of Multi-component Vesicle Membranes 
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In this work, we study the adhesion of multi-component vesicle membrane to both flat and curved 
substrates, based on the conventional Helfrich bending energy for multi-component vesicles and 
adhesion potentials of different forms. A phase fleld formulation is used to describe the different 
components of the vesicle. For the axisymmetric case, a number of representative equilibrium vesicle 
shapes are computed and some energy diagrams are presented which reveal the dependence of the 
calculated shapes and solution branches on various parameters including both bending moduli and 
spontaneous curvatures as well as the adhesion potential constants. Our computation also confirms 
a recent experimental observation that the adhesion effect may promote phase separation in two- 
component vesicle membranes. 



I. INTRODUCTION 

Adhesion is a fundamental step for many biological 
processes such as exocytosis, endocytosis. Cell adhesion 
also plays important roles in drug designs and drug deliv- 
eries as well as many biosensor applications P, 0] . There 
have been many experimental and theoretical studies fo- 
cusing on this subject 0-0] ■ While many of the past 
studies on the vesicle-substrate adhesion have focused on 
the case of a flat substrate [§-[1^. there have also been 
some works that address the complexity of curved sub- 
strate. For instance, theoretical and experimental stud- 
ies on the binding of a vesicle membrane to micro or 
nano-particles, or colloids have been conducted in [TBl - 
[igj . where the characteristic spherical substrates have 
radii much smaller than that of the vesicles. In the 
adhesion of a three dimensional vesicle to curved sub- 
strates has been studied where the curvature of the sub- 
strates are comparable to the curvature of the vesicles. A 
phase diagram for bound-unbound transitions has been 
presented. In this work, we study the adhesion of multi- 
component vesicle membranes to both flat and curved 
substrates. This is motivated by experimental studies 
of the modeled subjects. For instance in a recent experi- 
ment conducted by Gordon et al. it was observed that 
a mixed-lipid membrane can go through a local phase 
separation above critical dcmixing temperature due to 
its close proximity to a biological or non-biological sur- 
face. That is, adhesion can promote the phase separation 
for the mixed-lipid cell or vesicle membranes. 



In this work, we develop a phase field model to study 
the adhesion of multi-component vesicle membranes with 
a substrate through a specified adhesion potential. Fol- 
lowing our recent approach described in [20[ , we take the 
adhesion potential to be a function of distance between 
the membrane and the substrate. The strength of ad- 
hesion potential is considered to be distinct for different 
components. By minimizing the total energy of the sys- 
tem that includes bending energy, interfacial line tension 
and the adhesion energy, the equilibrium vesicle shapes 
can be computed for a variety of parameter values. As the 
initial attempt, we consider the case that both the vesicle 
membrane and the substrate are axisymmetric to sim- 
plify the computation. We present, in particular, a num- 
ber of typical equilibrium two-component axisymmetric 
vesicle profiles undergoing adhesion. The consistency be- 
tween the phase field description and its sharp interface 
limit is also briefly discussed. Moreover, a numerical ex- 
periment is conducted to support the conclusion of [J 
that the adhesion may promote phase separation for a 
multi-component membrane. 



II. MULTI-COMPONENT VESICLE 
MEMBRANE WITH ADHESION 

Equilibrium shapes of a multi-component vesicle are 
often described by minimizing an energy that includes 
elastic bending energy of the membrane and the line ten- 
sion energy at the interface between the components [2l[ . 
For the elastic bending energy of vesicle membrane, a 
common form adopted in the literature is that introduced 
by Helfrich [H: 
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where T is the membrane surface, H and K are the mean 
and Gaussian curvatures of T with A and b being the 
mean curvature bending modulus and the Gaussian cur- 
vature bending modulus respectively and a is the sponta- 
neous curvature. For simplicity, we consider the effect of 
mean curvature bending modulus and spontaneous cur- 
vature only. Thus b is set to be zero and Eq. [1] becomes 



Eb = J \{H 



In this paper, we focus on two-component vesicle mem- 
branes which have the liquid-ordered and the liquid- 
disordered phases, and the two phases have distinct bend- 
ing moduli and distinct spontaneous curvatures [23l - l25j . 
Let Fi and F2 be the parts of the surface representing the 
two different phases with Ai , A2 being their correspond- 
ing bending moduli and ai, 02 being their corresponding 
spontaneous curvatures, respectively. The total elastic 
bending energy for the two-component vesicle is 

Eb^ [ XiiH -aif dx+ [ A2(i/-a2)' dx. (2) 

The line tension energy, which is essentially an interfacial 
energy between the two phases is given as 



E, 



a dl , 



(3) 



where a is the constant line tension at the interface. Eqs 
dH) and ^ together define the total energy. 



E — Eh + El , 



(4) 



with a minimum of E describing the shape of an equilib- 
rium two-component closed membrane. 

Note that the vesicles or membranes discussed so far 
are free and not bounded to other objects. To incorporate 
the adhesive interaction with a substrate, an additional 
energetic contribution due to adhesion should be added 
to Eq. (gD: 

Etot.i= f Ai(ff-ai)2 dx+ / A2(ff-a2)2dx 

+ I (J dl - I PF(x) dx , (5) 



where 



X e Fi 
xe F2 



(6) 



is the adhesion potential which varies with respect to the 
position X on F = Fi U F2. In the above, wi and W2 are 
the corresponding strengths of the adhesion potential ex- 
perienced by the liquid-ordered and the liquid-disordered 
phases. A representative form of W is that of a Gaussian 
form given by. 



'u;icxp(-d(x)2/e2), 
W2 cxp(-d(x)2/e2), 



X e Fi 

xeF2 



(7) 



where d(x) is the distance from x to a flat/curved sub- 
strate, and e is a small number. Notice that when e 
approaches zero, the adhesion potential converges to a 
sharp contact potential, a scenario that has been investi- 
gated in earlier studies 0, Q . While we use the Gaussian 
potential ([7]) in most of this paper, to offer a comparison, 
we also consider the Leonard- Jones type potential, 
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which induces a narrow repulsive region between vesicles 
and the substrate. The constant P and the exponent a 
determine the thickness of the repulsive region and the 
rate of change of the adhesion potential, respectively. 



A. Phase field formulation 

To be able to effectively describe the different phases 
of the two-component vesicle, we use a phase field for- 
mulation which has become very popular in recent years 
in the modeling and simulations of vesicle deformations 
[20l . [26l - [3]| . A phase field function can be used to describe 
the vesicle with the phase field bending energy as formu- 
lated in [2^, [231 . Adhesion energy can be incorporated 
into the phase field formulation as shown in [20|. For 
multi-component vesicles, order parameters can be used 
to describe both the vesicle and its two components [s^] . 
On the other hand, for a vesicle with a fixed topology, 
one can also use a direct (explicit) surface representa- 
tion for the vesicle along with an order parameter (phase 
field function) to describe the two different phases of the 
membrane (so. Issj . For the axisymmetric case considered 
here, it is particularly effective to adopt a sharp interface 
representation of the vesicle surface given by the revolu- 
tion of a simple one-dimensional curve with an arc-length 
parametrization and a phase field representation of the 
different phases on the vesicle which is also a function of 
the arc-length. 

Specifically, let F be the vesicle surface, a phase field 
function 77 = 77(x) is introduced over F which may be 
used to represent either a material composition profile 
or a fictitious density of the lipids on the surface of the 
membrane and distinguishes between the liquid-ordered 
and liquid-disordered phases. As an illustration, we focus 
on the latter case so that in the liquid-ordered phase, 77 is 
specified to be and is colored as blue in figurc[TJ in the 
liquid-disordered phase, 77 is assigned to be —1 and col- 
ored as red. In the interface between the liquid-ordered 
and liquid-disordered phases, rj rapidly, but continuously, 
changes from +1 to -1. Note that the phase field regular- 
izes the sharp interface between the two different phases 
into a diffused one, and thus provides a more general de- 
piction of the two-component vesicle in both the mixed 
and de-mixed states. The total energy ([5]) of the model 
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in terms of the phase function rj is given by 



||Vr?7|' + ^(??) 



dx 



w(co + C277)^(d(x)) dx, (9) 



where the first and third terms are phase field formulas 
for elastic bending energy and adhesion energy, respec- 
tively. The term cq + cirj is a phase field representation 
of Ai and A2 . When x is away from the interfacial region, 
co + ci = Ai and cq — ci = A2. The liquid-ordered phase is 
stiffer than the liquid-disordered phase, hence we always 
assume ci > 0. The term a(co + c^rf) is considered as the 
phase field analog of spontaneous curvature, and 

a(co + C3) = ai, a(co - C3) = 02 

if X is away from the interfacial region. Similarly, 
w(cq + C2?7)P(d(x)) is viewed as an approximation of the 
adhesion potential VF(x), with 

w(co -f C2) = Wl, w(co - C2) = W2 

when X is far away from the interfacial region. The sec- 
ond term is a phase field approximation for the line ten- 
sion energy where a double well potential function 



(10) 



is incorporated. Vr?? is the surface gradient of 77 which is 
the projection of V77 onto the tangent plane of V. Notice 
that to make Vr?? well defined, the function 77 should be 
defined away from the membrane such that d?7/dn ~ 
where n is the normal vector of T. The enclosed vol- 
ume and total area of the membrane are assumed to be 
invariant. Meanwhile, the total amount of lipids is con- 
served. Thus three constraints are imposed during the 
minimization of the total energy ([5]): 

y" dx = A, y" = Vol, j 77(x) dx = C. (11) 

The constraint /p 77 (x) dx = C denotes the difference in 
the surface areas of the two phases in the sharp interface 
limit. 



B. Axisymmetric Setting 

In the present work, we focus on the axisymmetric 
membrane adhered on a flat/curved substrate. In this 
setting, the membrane surface is determined by evolving 
a 2-d curve. A vesicle with a flat substrate is schemati- 
cally shown in figure [T] The axisymmetric vesicle surface 
is generated by evolving a curve parameterized by the 
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FIG. 1. Schematic diagrams of axisymmetric two-component 
adhered vesicle membranes. Blue and red colors indicate the 
liquid-ordered and the liquid-disordered phases, respectively, 
s = So specifies the phase boundary. Away from the interface, 
the blue phase has mean curvature bending modulus Co -I- ci 
and adhesion potential w{co + C2)P{d{x)), while the red phase 
has mean curvature bending modulus Co — ci and adhesion 
potential w(co — C2)P(d(x)). 

arc- length s, and the total length of the generating curve 
is denoted by s. The flat substrate is located at z = 
with the two different phases being distinguished by the 
red and blue colors. The transition point from red to 
blue is located at s = sq- The figures on the left and right 
show different configurations with the blue phase and the 
red phase being adjacent to the substrate, respectively. 
For easy reference, we denote the left as red-blue vesicle 
membrane, and the right as blue-red vesicle membrane. 
The tangent angle (j) is measured from the radial direc- 
tion and r is the distance of a point on the membrane 
from the axis of symmetry. 

The mean curvature of the vesicle can be explicitly 
expressed by r and cj) as 



1 / , sin ( 



where prime represents the derivative with respect to arc- 
length s. The line tension energy term in ^ becomes 

2na + j-^iv' - 1)' ds. 

We nondimensionalize all the parameters and choose 
Co to be 1. Then the phase field model ^ is reduced to 



E{r]) = 27r / (1 + cirj) [H - a{l + C377)] ' r ds 
Jo 



-27rcr 



ds 



2 ' 4^ 

-27r / u;(l + C2?7)P(d(x)) r ds. (12) 
Jo 

Additionally, the constraints (jlip in the axisymmetric 
case lead to 



(cosT)' 



(13) 
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with T being the arc-length parameter of the reference 
unit sphere, 



and 



TT / r'^z'ds = Vol, 
/o 



ri[s)rds = C . 



(14) 



(15) 



The pointwise constraint (jl3p . which is in fact equivalent 
to the one in (jlip . is referred as the lateral incompress- 
ibility condition for the membrane [s^ ]. 

The shape of the membrane is determined by minimiz- 
ing the total energy p^ . subject to constraints (fT3)l . ((M]) 
and (fT5|) . It satisfies the Euler-Lagrange equations given 
by: 



[H" 



r'H' 



2H{H^ - K) + 2aHH{\ + cgr/)] - 2{^iH + p + t7]H) 



2aH 



dP (5(d(x)) 
d(d(x)) (5n 



2HPj = (16) 



ci [F - a(l + C3?7)] -2ac3ff + + 

I 



r 7/ 
r 



T — C2WP = , 



(17) 



where H = {1 + cirj) [H — a(l -I- csry)] and and r are 
the three Lagrange multipliers for the three constraints, 
respectively. The other equations from geometry are (25l . 
Hi: 



4)' = 2H - 



sm( 



r = cos ( 



sm( 



Boundary conditions are imposed as follows: 

H'{Q) = H'{s) = 0; r(0) = r(s) = 0; 
0(0) = 0, <j){s) = tt; 77' (0) = 77' (J) = 0. 



III. NUMERICAL EXPERIMENTS 



(18) 



(19) 



We numerically solve Eqs (|T6|) through ((18)1 subject to 
boundary condition (HI]) using the MATLAB ODE solver 
BVP4C. We tested the convergence of the computed re- 
sults in our numerical simulation. 



A. Adhesion with Gaussian potential 

We first present the numerical results of a few adhered 
vesicles using the Gaussian adhesion potential. 

In Figure [21 several numerical solutions for Eqs. (fTBl) 
through (|19p are shown. The curves are the cross-sections 
of the vesicle shapes with the blue-colored region repre- 
senting the liquid-ordered phase and the red-colored re- 
gion representing the liquid-disordered phase. The cor- 
responding parameter values used in the experiments 
are taken as area=47r, volume=3.5, C = —0.3, a = 2, 
ci = 0.4, C3 = 0, a 0, C = 0.01, and e = 0.15. Note 




(B) 





FIG. 2. Membrane shapes for different w. Area=47r, vol- 
ume=3.5, area difference C = —0.3, a = 2, ci = 0.4, C3 = 0, 
a = 0, ^ = 0.01, e = 0.15; (A) C2 = -0.05, w = 1.0, 1.5, 2.0, 
from left to right; (B) C2 = 0.2, w = 2.0, 3.0, 5.0, from left to 
right. 

that ci — 0.4 implies that the ratio A1/A2 ~ Abiue/Arod of 
the mean curvature bending moduli between the blue and 
red phases is 1.4/0.6. Moreover, a = and C3 = imply 
that both phases have zero spontaneous curvature. Blue- 
red adhered vesicles on flat substrate are shown in figure 
[2]-A where C2 = —0.05, that is, Wi/w2 = Whiue/wred in 
([7]) is equal to 0.95/1.05. Except the adhesion-associated 
parameter w, all the parameters are kept fixed as speci- 
fied. Some adhered shapes of red-blue vesicles are shown 
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FIG. 3. Convergence of adhered membranes as e 0. 
e = 0.15,0.10,0.05,0.01, from left to right. Area=47r, vol- 
ume=3.5, area difference C = —0.3, ct = 2, ci = 0.4, 
C2 = -0.05, C3 = 0, a = 0, C = 0.01, w ^ 2. 

in figure 2-B for various w and C2 = 0.6. Notice that 
there is a shght protrusion of vesicles into the flat sub- 
strate. This is due to the lack of repulsive effect in the 
Gaussian form of the adhesion potential. 

Convergence of the adhered vesicle shapes as e ap- 
proaches zero is presented in figure [S] Theoretically, by 
the standard asymptotic analysis, rj tends to converge to 
tanh ( ) when ^ approaches zero, where sq indicates 
the location of the interface which depends on the area 
difference C (see appendix A). Numerically, for relatively 
larger e, the shapes of two-component vesicle membranes 
protrude into the flat substrate more significantly. As 
e gets smaller, the protrusion becomes less visible and 
finally vanishes in the limit of e — 0, that is the case 
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FIG. 4. Energy comparison between blue-red vesicle and red- 
blue one. w = 2, ci — 0.3 are fixed. C2 varies from -0.05 to 
0.6. The solid curve represents energy versus C2 for the red- 
blue vesicle, while the dash curve is for the blue-red vesicle. 
The shapes of the vesicles at the four end-points of the curves 
are shown. Area=47r, volume=3.5, area difference C — —0.3, 
cr = 2, C3 = 0, a = 0, ,t = 0.01, e = 0.15. 
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FIG. 5. Transition curves for w=1.5, 2.0, 3.0. Above each 
curve, red-blue membrane is more stable; the blue-red mem- 
brane is more stable below the curves. The other parameter 
values are: Area=47r, volume=3.5, area difference C = —0.3, 
a = 2, C3 = 0, a = 0, ,t = 0.01, e = 0.15. 

corresponding to the effective contact potential [1, 0, [l^l 

W{x.) I Q otherwise. ^"^^^ 

Figure 2] shows the energy comparison between blue- 
red vesicle membrane and red-blue one under the influ- 
ence of adhesion. Here area=47r, volume=3.5, C = —0.3, 
CT = 2, ci = 0.3, C3 = 0, a = 0, ^ = 0.01, e = 0.15, 
w = 2.0 are fixed. The two shapes on the left correspond 
to parameter C2 = —0.05 with wuue/wied = 0.95/1.05. 
The blue-red vesicle, which is more deformed from the 
free vesicle shape, has lower energy and is more stable 
than the red-blue one in this case. On the other hand, 
the two shapes on the right correspond to C2 = 0.6 lead- 
ing to Wbiuc/w^rcd = 1.6/0.4. There, the red-blue vesicle 
which is deformed more from the free shape has lower en- 
ergy and is more stable than the blue-red one. In general, 
for C2 < 0, uifod is larger than Wbiuc and blue- red mem- 
brane is more deformed from a free shape. Whereas, for 
C2 > 0, Wbiuo is larger than z/;,od and red-blue membrane 
is more deformed from a free shape. 

The energy comparison of the four shapes as discussed 
above may lead us to believe that the vesicle whose com- 
ponent adjacent to the substrate endures stronger adhe- 
sion is more stable. However, this is not always the case. 
As seen in figure IH there is an anomalous region with 
C2 between zero and 0.1217, and the blue phases suffer 
stronger adhesion but the blue-red vesicle is more stable. 
Outside this region, the shape with stronger adhesion on 
the component adjacent to the substrate is more stable. 
We observe that the existence of the region (0, 0.1217) 
is due to the nonzero values of ci and C which model, in 
the two-component system, the effects due to the differ- 
ences in the bending moduli and the surface areas of the 
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FIG. 6. Dependence of C2 on w for phase transition when 
ci — 0.3 and 0.0, respectively. Area=47r, volume=3.5, area 
difference C = -0.3, a = 2, C3 = 0, a = 0, C = 0.01, e = 0.15. 

two phases. When both ci and C approach zero which 
is the hmit where both components have the same bend- 
ing moduh and equal surface areas, the multi-component 
vesicle reduces into a single component vesicle. And the 
anomalous region shrinks to the point C2 = 0, also the 
transition point converges to C2 = 0, 

In addition, the transition point in figure 2] strongly 
depends on ci and the adhesion potential w. The depen- 
dence of the transition points on ci is shown in figure [S] 
For various w — 3.00,2.00,1.61, we plot the transition 
curve ci versus C2. The red-blue vesicle, corresponding 
to the parameter pair (ci,C2) above each curve, is more 
stable; while the blue-red vesicle is more stable if the 
parameters ci , C2 are chosen from the region below the 
transition curves. 

Figure [5] shows that the transition always occurs for 
C2 > 0. One may wonder if this is universally true for 
any adhesion potential w. The answer is provided via 
figure [6l the transition curve w v.s. C2 for fixed ci. In 
figureO the dash curve corresponds to ci = 0.3, and the 
solid curve corresponds to ci =0. Then we obtain the 
dependence of C2 on it; when the transition occurs. For 
ci = 0.3, w varies from 1.61 to 3; while for ci = 0.0, w 
varies from 1.26 to 3. Above the transition curve, the red- 
blue vesicle is the more stable one; while below the curve, 
the blue-red vesicle is more stable. The bound blue-red 
vesicle with given parameters will break away from the 
adhered state around w = 1.60(1.25) when ci = 0.3(0.0) 
and change to an unbound (free) vesicle. We can thus 
claim that C2 is always positive for any bound adhesion 
potential w when ci is fixed. 

We also take the effects of the spontaneous curvature 
and the substrate curvature into consideration [3^, [IJ ■ 
In figure[7]thc fiat substrate is now replaced by a concave- 
up spherical substrate with radius R = 3. Here area=47r, 
volumc=:2.7, area difference C = 0, cr = 1, ci = 0, C2 = 0, 
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E=8.6450 



FIG. 7. Effect of spontaneous curvatures on the vesicle 
shapes. Representative shapes of vesicle membranes are 
shown. Area=47r, volume=2.7, area difference C = 0, a — 1, 
ci = 0.0, C2 = C3 = 0.5, a = 2/3, C = 0.01, e = 0.05, w = 2. 

C3 = 0.5, a = 2/3, £, = 0.01 and e = 0.05. We take w = 
and w = 2 respectively to find both free vesicles and ad- 
hered vesicles. Differing from the previous figures in this 
subsection, the red and blue colors in this experiment 
indicate the phases having different spontaneous curva- 
tures. The blue phase has a bigger spontaneous curva- 
ture a(l -|- C3) = 1 while the red phase has a smaller one 
a(l — C3) = 1/3. Al and Bl show the free vesicles with 
specified parameters. Notice that the red phase, which 
possesses a smaller spontaneous curvature, has shapes 



0.08 
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FIG. 8. Dependence of C2 on C3 for phase transition. Rep- 
resentative shapes of vesicle membranes are shown above or 
below the transition curve. area=47r, volume=2.7, area dif- 
ference C = 0, (7 = 1, ci = 0.0, a = 2/3, C = 0.01, e = 0.15, 
It) = 2. 
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less elongated than the blue phase. A2 and B2 are ad- 
hered blue-red and red-blue vesicles. By comparing the 
total energy of these two adhered vesicle, we find out that 
the blue-red vesicle with the red phase at bottom, whose 
spontaneous curvature matches with the curvature of the 
curved substrate, is more stable than the red-blue vesicle. 

We choose here a volume with value 2.7 which is 
smaller than the previously chosen value of 3.5 because 
larger osmotic pressure difference makes the effect of 
spontaneous curvatures on the vesicle shapes more visi- 
ble. 

Similar to the discussion of stability transition in fig- 
ure [5] and figure [6l when the spontaneous curvature effect 
is taken into account, we can also find the dependence 
of transition points on ci, C2, C3, w. Without repeatedly 
showing many figures, we present here only the transition 
curve C3 v.s. C2 in figurc[8]with ci = 0, w = 2 fixed. Other 
parameters are set as area=47r, volumc=2.7, area differ- 
ence C = 0, a = 1, a = 2/3, ^ = O.Of, e = 0.f5. Above 
the transition curve, a representative red-blue vesicle is 
shown which is more stable than the blue-red vesicle; 
while below the transition curve, the blue-red vesicle is 
more stable. 





FIG. 9. Membrane shapes for Leonard- Jones type potential. 
(A) w = 2, /? = 0.05, 0.035, 0.025, from left to right; (B) f3 = 
0.05, w = 1.2, 2.0, 2.8, from left to right. Other parameters: 
Area=47r, volume=3.5, area difference C = —0.3, cr = 2, ci = 
0.3, C2 = -0.05, C3 = 0, a = 0, C = 0.01, a = 5 



B. Adhesion with Leonard- Jones potential 

When the Gaussian adhesion potential is used, as seen 
in the above, slight protrusion of vesicles into the sub- 
strate may occur. If this protrusion is of any practical 
concern, one possible way to avoid the protrusion is to 
replace the Gaussian adhesion potential by a Leonard- 
Jones type potential. We demonstrate, by employing the 
following potential. 



P(d(x)) = -4 



d(x) 



(3 



d(x) 



a/2-i 



(21) 



that the shapes of adhered vesicles without any protru- 
sion can be obtained using our phase field formulation. 
The key difference between the Gaussian type potential 
([7]) and the Leonard- Jones type potential (|21|) is that 
Gaussian potential is globally attractive, while there is 
a narrow region with d(x) between zero and [3 where 
the Leonard- Jones type potential is repulsive. Such a 
repulsive region can prevent the vesicle membrane from 
protruding into the substrate. 

An axisymmetric simulation is shown in figure [H] where 
we take a = 5. In figure [HI^A) , w = 2 is fixed, and the 
parameter (3 is decreased from left to right. We can see 
that the repulsive region between substrate and vesicle 
membrane narrows down. The shape deformation under 
the influence of adhesion is shown in figure [5KB) where 
the repulsive region is fixed but w is increased from left 
to right. The stability transition curves, similar to that 
obtained previously for Gaussian potential, can also be 
obtained for Leonard- Jones type potential. Given the 
similarities in the findings, we do not repeat the discus- 
sion here. 



C. Promotion of phase separation 

In this subsection, two numerical experiments are pre- 
sented to support Gordon and coworkers' experimental 
observation that adhesion may promote the phase sepa- 
ration The experiments involve the phase field sim- 
ulations for relatively large diffuse interface parameter 
^ which is consistent to the experimental setting. The 
Leonard-Jones type adhesion potential is employed in 
this subsection. 

In the first experiment, the function 77(3) is viewed as 
a chemical composition function, which can be consid- 
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FIG. 10. Adhesion induces phase separation. An almost ho- 
mogeneous free vesicle is phase-separated by a Leonard- Jones 
adhesion. 
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ered as a compositional fluctuation around a homoge- 
neous state with composition 779. The constraint (jl5p is 
specified as 

/ ri{s)rds = / rj^rds (22) 
Jo Jo 

We now demonstrate that for an equilibrium free vesi- 
cle with associated 77(5) almost homogeneous, namely 
r]{s) « Tjo, after adding the adhesion, it displays the phase 
separation behavior with r][s) changing into a tanh-like 
profile representing two distinct phases. 

The numerical experiment is presented in figure 1101 
In (Al-2), for area = 47r, volume = 3.5, 779 = — l/\/3,C = 

0.33, cr = 5/(1-0.8770), Ci = 0.8/(1-0.8770), C2 = 1.0/(1- 

I.O770), C3 = 0, a free vesicle is computed first. The chem- 
ical composition function 77 and the associated vesicle 
profile shown in (Al-2) correspond to the stable energy 
minimum (in the absence of adhesion). The change in 
composition is relatively small, representing a state of 
mixed phases in much of the vesicle. By adding Leonard- 
Jones adhesion with w = 5(1 — 1.0r/o)/(l — 0.8770),/? = 
0.023, a = 5, an adhered vesicle with associated 7] corre- 
sponding to the new energy minimizcr is shown in (Bl- 
2). A phase separation occurs with a much more sig- 
nificant phase difference max(77) — min(7/) = 1.9703 in 
(B2). In comparison, we have max(77) — min(77) = 0.2345 
in (A2). which is only 11.90% of the phase difference 
for adhered vesicle. We thus see that adhesion can sig- 
nificantly promote phase separation. Furthermore, the 
bending modulus of the phase-mixed vesicle Amixod is 
roughly around 1/(1 — O.8770), and the bending moduli of 
the phase-separated vesicle are Abiuc = 2.2649/ (1 — 0. 8770) 
and Arod = 0.6873/(1 — O.8770). These parameter values 
imply that Arod ~ Amixod and Abiue/Arod = 3.2954 which 



(A1):w=0 (B1):w=2 (C1):w=10 




r r r 




0.5 10 0.5 10 0.5 1 



TIT 

FIG. 11. Promotion of phase separation. As w gets bigger, the 
two phases of 77 near the interfacial region are separated more. 
Area=47r, volume=3.5, area difference C = 1.7, ^ = 0.13, 
a = 1, ci = 0.05, C2 = -1.00, C3 = 0, C = 1, /3 = 0.05, a = 5. 



agree with the ratios of the bending stiffnesses of the or- 
dered, the disordered and the mixed phases available in 
the literature 0, [13, [H, [H] . 

The effect of adhesion on phase separation can be fur- 
ther demonstrated in the next experiment. Here again, 
we take the order parameter 77(5) to be a labeling func- 
tion for relative lipid density which stays within [—1, 1]. 
To model a wider diffuse interfacial layer corresponding 
to a relatively large ^ and yet maintain the bound on rj, 
we replace the double well potential function (1 101) by the 
following double obstacle potential function [39j : 

(l-f77)ln(l+7/) + (l-77)ln(l-77)4^7(l-772)-21n2 . (23) 

where 7 = 1-1-2 In 2 is a constant describing the height 
of the potential barrier. 

Figure [11] shows a numerical experiment while we 
choose ^ = 0.13 so that for unbounded vesicles, there 
would be a wide interfacial regions with less dramatic 
phase separation effect. Then for various w = 0, 2, 10, we 
compute the corresponding adhered vesicle shapes and 
the corresponding lipid density functions 77. We high- 
light the densities near the interfacial region (the boxed 
portion of the vesicle profile). With the plots using the 
same scaling in the arc-length, one can see that with a 
larger adhesion strength w, the interfacial layer gets nar- 
rower and the separation between the red and blue phases 
becomes sharper and more dramatic. 

IV. CONCLUSION 

In this work, we develop a phase field model for the 
adhesion of the multi-component vesicle membrane to a 
flat/curved substrate. Some representative vesicle shapes 
are presented. The infiuence of C2, which measures the 
contrast of the adhesion effect for the different phases, 
on the stability of membrane shapes is discussed here. 
It turns out that the multi-component vesicle, whose 
lower part suffers stronger adhesion, is more stable in 
most cases. We also consider the the phase transition 
from blue-red vesicle to red-blue vesicle, and the influ- 
ence of other parameters such as the relative contribu- 
tion of the adhesion and the bending energy. The effect 
of spontaneous curvature is also numerically observed by 
examining the adhesion of the vesicles on a curved sub- 
strate. We also present vesicle shapes corresponding to 
the Leonard-Jones type adhesion potential to show that 
the repulsive region can effectively prevent the protru- 
sion of the vesicles into the substrates should this be of 
practical concern. Finally, we numerically examine the 
fact that adhesion can promote the phase separation for 
multi-component vesicle membrane. 

Although the numerical studies presented here are fo- 
cusing on the axisymmetric configuration, the phase field 
formulation is applicable to more general settings such 
as in the study of the interaction between the multi- 
component vesicle and a patterned substrate. It can also 
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be used to extend the study of membrane-mediated parti- 
cle interactions [l8l.[l9| to multi-component vesicles. The 
present phase field formulation of the multi-component 
vesicle and substrate interaction can also be extended 
in several directions. For instance, by incorporating the 
Gaussian curvature contributions to the bending energy, 
fusion of multi-component vesicles can be studied. By 
adding an entropic contribution to the free energy, we 
can also consider the weak adhesion regime where fluctu- 
ation of the vesicle shape due to thermal excitation plays 
an important role. 



Appendix A: Asymptotic analysis for 77 

In section Hill the convergent behavior of the adhered 
vesicles as the interfacial width e approaches zero is 
briefly mentioned. It is claimed that the phase field 
function 77 approaches tanh(^^) as the parameter ^ 

approaches zero. A boundary layer calculation is carried 
out to support such an observation [40| . 

The lipid density function r](s) is governed by the Eq 
(jl7|l and can be rewritten as 

C(ci [H -a{l + C3T])] ^ - 2ac3H + t - C2wP 



(Al) 



An asymptotic analysis in the sharp interface limit is 
carried out as follows. For simplicity, we only consider 
the 0(1) terms for outer and inner layers here. The in- 
ner layer and the outer layer are the regions around the 
interface and away from the interface, respectively. 

Within the outer layer, rj does not change much with 
respect to the arc-length and we expand 77(5) as 



vis) ~ ^(s) + ??7i(s) -f 



(A2) 



We substitute (|A2[) into (jAip and obtain an equation for 
?7o at the lowest order: 



The solutions are 



(77o'-l)%r = 0. 



770 = 0, ±1 



and we choose 77 = 1 and 77 = —1 in the two sides of the 
interface to represent the blue and red phases, respec- 
tively. Here let us assume t] = —1 in the side s < sq; and 
77 = 1 in the side s > sq. 

In the inner layer, we expect 77 to vary rapidly from —1 
to 1 and introduce a new stretched arc-length variable 



So 



where fx is an yet undetermined parameter. We now ex- 
press r; as a function of the stretched variable S and ex- 
pand ?7 as 

?l{S)^?,oiS)+^miS) + -- - . (A3) 

Notice that 

r{s) = r{so + eS) = r{so) + ^Sr'iso) + 0{en, 



Upon rewriting Eq. (|Aip in terms of the stretched vari- 
able S and, subsequently, using Eq. (jA3p we obtain 

+ - l)?7o[r(.so) + eSr'is,)] + 0{e) = (A4) 

There are two possible choices for {i. If /i balances the 
second and third terms, namely, 2 — /i ~ 0. Then the 
leading order term is: 



-riso) 







which implies 770 = aS + b. However, this solution does 
not satisfy the matching conditions 

77o(-oo) = 770(0) = -1, 77o(+oo) = 77o(s) = 1 . 

Another choice is to balance the first and third terms 
of equation (|A4p . namely, yu = 1, and one gets the leading 
order term: 

-''(so)^^ + {V^- l)mr{sa) = 0. 

If the boundary condition 7fo(— 00) = —l,rjo{+co) ~ +1 
are imposed, a typical solution of the above nonlinear 
equation, which satisfies the matching condition, is 



770(5) = tanh 



(9^ 



Then the composite solution of (jAip . given by 
inner solution + outer solution - matching solution, 

takes the form 

^o(5) + 77o-(s) - 7,o-(0) + 77+ (.s) - 77+ (s), 

and explicitly 

' s - So ^ 



77(5) ~ tanh ^- 
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